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ABSTRACT 



We have performed an On-The-Fly (OTF) mapping survey of 
^^CO(J = 1 — 0) emission in 28 Virgo cluster spiral galaxies using the Five 
College Radio Astronomy Observatory (FCRAO) 14-m telescope. This survey 
aims to characterize the CO distribution, kinematics, and luminosity of a large 
sample of galaxies covering the full extents of stellar disks, rather than sampling 
only the inner disks or the major axis as was done by many previous single 
dish and interferometric CO surveys. CO emission is detected in 20 galaxies 
among the 28 Virgo spirals observed. An atlas consisting of global measures, 
radial measures, and maps, is presented for each detected galaxy. A note 
summarizing the CO data is also presented along with relevant information 
from the literature. The CO properties derived from our OTF observations are 
presented and compared with the results from the FCRAO Extragalactic CO 
Survey by Young et al. (1995) which utilized position-switching observations 
along the major axis and a model fitting method. We find that our OTF derived 
CO properties agree well with the Young et al. results in many cases, but the 
Young et al. measurements are larger by a factor of 1.4 - 2.4 for seven (out of 
18) cases. We will explore further the possible causes for the discrepancy in the 
analysis paper currently under preparation. 

Subject headings: atlases — galaxies: clusters: individual (Virgo) — galaxies: 
ISM — radio lines: galaxies 
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Introduction 



The high galaxy density and the proximity make the Virgo cluster a particularly inter- 
esting laboratory for a galaxy evolution study. Its dynamical evolutio n is still in progress , 



and evidence for significant environmental effects is ubiquitous (e.g., iChung et al.l 120071 ). 
The nearness of the Virgo cluster makes it possible to observe its member galaxies with 
excellent spatial resol ution (1" ^ 90 pc). It is the first cluster for which significant 



HI imag ing was done (Ivan Gorkom et al.l Il984f). and v arious new HI surveys such as AL- 
FALFA JciovaneUi et allboOTh and VIVA f lchunelboOTl l have recently been conducted. The 
Virgo cluster has also been the subject of many recent mu lti- wavelength surve ys, such 
as in radio continuurn by NRAO VLA Sky Sur vey (NVSS; Condon et al. IQQsh. in Ha 
jKoopmann fc KennevI [20oi : Ichemin et al.l[2006[ ) , in UV by FAUST JSrosch et al.l Il997h 



and Galaxy Eyoluti on Explorer (GALEX) (e.g., iDale et al.l l2007l ). and in IR by Spitzer 



( Kenney et al. 20081 ). High quality optical ir nages are also ay ailable from the Hubble Space 



Telescope (HST)/ACS Virgo Cluster Survey (ICote et al.ll2004l ) and Sloan Digital Sky Survey 
(SDSS). 

Here, we present the results from a new imaging survey of J = 1 -^^ ^^CO emission 
in a large sample of Virgo galaxies in order to address the distribution and characteristics 
of dense molecular gas in these galaxie s . It is well established that molecular clouds are the 
sites of ongoing star formation (iLarsonl (120031 ) and references therein), and carbon monoxide 
(CO) is the most commonly used tracer of molecular hydrogei i (H2), which is the mos t 



abundant but invisible component of cold and dense clouds (e.g.. iSolomon &: Barrettlll99ll ). 
The global H2 content and its distribution in galaxies and a comparison with other gas and 
stellar components as a function of morphological type, luminosi ty, and environment are 
some of the key insights one can derive from CO observations (see lYoung &: Scovilld (Il99ll ) 
and references therein). 

Existing CO surveys of Virgo cluster galaxies suffe r from limited spat ial coverage and 
small sample sizes. The FCRAO extragalactic survey (lYoung et al.l Il995l ) is the first CO 
survey covering a wide range of distance, diameter, morphological type, and blue luminosity 
of some 300 galaxies conducted using the FCRAO 14-m telescope, including CO measure- 
ments (detections and upper limits) of 65 Virgo galaxies. Because of the large observing time 
required, these observations were conducted in the position-switching mode, primarily along 
the optical major axis of the disks, and the global CO line luminosity was derived assuming 
a model distribution. More recently, h igh angular resolution CO i mages have been o btained 
using interferometric measurements by lSakamoto et al.l (119991 ) and lSofue et al.l (120031 ). These 
measurements reveal a detailed molecular gas distribution at 100 pc scales, but they are lim- 
ited only to the central region of galaxies. The BIMA SONG (Survey Of Nearby Galaxies; 
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Heifer et al.l (120031 )) has also carried out an imaging survey of CO emission in several Virgo 
spiral galaxies, incorporating the short spacing data from the NRAO 12-m telescope. The 
total number of Virgo galaxies imaged by the BIMA SONG is small (six), however. 

These earlier CO observations have revealed important insights on the molecular in- 
terstellar medium (ISM) in these galaxies. For example, CO emission is often centrally 
concentrated, in contrast to the centrally deficient HI distributions commonly found in these 
galaxies. The molecular gas dis tribution also shows lit tle evidence for any infl uence of the 
gas s tripping mechanisms (e.g., Kennev fc Younj 19891). No central CO peak ( Heifer et al. 



20031 ) or nuclear molecular rings (e.g. 



lono et al 



20051 ) are seen in other cases. The in- 
fluence of cluster envir onment on the molecular content is still poorly understood (e.g., 
BoseUi fc Gavazzilbooeh . 



A distinguishing characteristic of our new CO survey is the complete imaging of ^^CO 
(J=l-0) emission of a large sample (28 galaxies) of Virgo spirals using the On-The-Fly (OTP) 
mapping mode of the Five College Radio Astronomy Observatory (FCRAO) 14-m telescope. 
Our map size of 10' x 10' is larger than the optical diameter D25 and it covers the entire 
stellar disk of each galaxy. The 45" angular resolution of the new CO images is well matched 
to the existing VLA HI data and is well suited for comparison with other high resolution 
multi-wavelength data. The specific questions we aim to address are : 

1) To what extent is the CO distribution governed by the disk dynamics? 

2) Is there a clear phase transition between HI and H2 as a function of the interstellar 
radiation field? 

3) Are there any systematic differences in the CO properties of spiral galaxies in different 
environments? 

We present in this paper the data and the CO atlas from our OTF mapping survey. 
In section 2, we describe the sample selection. Observations and data reduction process are 
described in Section 3. The CO atlas and CO prop e rties are presented in Section 4, and 
our results are compared with those of lYoung et al.l (119951 ). Molecular gas distribution in 
individual galaxies is discussed in Section 5, and the summary and conclusion are given in 
Section 6. The analysis and interpretation of the data addressing the above questions will 
be presented in our subsequent papers. 
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2. THE SAMPLE 

We initially selec t ed th e 42 Virgo spiral galaxies in the magnitude limited sample studied 



by iKenney fc Youngl (Il988l ). The primary goal of our project is investigating the spatially 
resolved distribution of molecular gas and their relation to other tracers of activities among 
the Virgo spirals, and a large angular size is an important consideration. We also considered 
the total CO flux to increase the detection rate. Given the observing time limitations, a 
subset of 28 galaxies was observed in the order of CO strength (Sco > 200 Jy km s^^) and 
included in this study. The observed galaxies have an optical diameter of 3-10 arcminutes 
and span a wide range of morphology (Sa to Sc), surface brightness (9 < BSp < 12), and dust 
mass (10^ < Mdust ^ IO'^Mq). The basic properties of the observed galaxies are summarized 
in Table 1 and are illustrated in Figure [H 



^The E-SO galaxy NGC 4649 is included in our sample because it is a companion of a late type spiral 
NGC 4647. 



Table 1. General Properties of the Sample 
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Table 1 — Continued 
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Note. — The columns are as follows. 
Col. (1) — NGC number 

Col. (2) & (3) — Right Ascension and Declination at epoch J2000 



Col. (4) — Morphological Type (RSA : Revised Shapley-Ames Catalog) 

Col. (5) — Numerical Morphological Type code from LEDA (Lyon/Meudon Extragalactic Database 
http: / /leda.univ-lyonl .fr) 

Col. (6) — Angular distance from M87 {a = 12h30m49.18s, 6 = 12d23m28.9s) 
Col. (7) — Major axis diameter at 25th mag arcsec"^ in the 5-band (LEDA) 

Col. (8) — Total S-band magnitude corrected for galactic extinction, internal extinction, and k-correction 

(LEDA) 

Col. (9) — Dust temperature from Young et al. (1989) 
Col. (10) — Inclination (LEDA) 

Col. (11) — Mean heliocentric radial velocity (LEDA) 



I 
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3. OBSERVATIONS & DATA REDUCTION 
3.1. Observations 

We carried out On-Tlie-Fly (OTF) mapping observations of CO emission in 28 Virgo 
galaxies over several observing sessions between January 2002 and February 2003 using 
the SEQUOIA (SEcond QUabbin Optical Imaging Array) focal plane array receiver on the 
FCRAO 14-m telescope. The SEQUOIA consists of 16 horns, each with 45" beam size, 
configured in a 4x4 array. The backend system used is the Quabbin Extragalactic Filterbanks 
(QEF), which consists of sixteen independent spectrometers each with 64 channels at 5 MHz 
resolution, resulting in a total bandwidth of 320 MHz {AV ~ 830 km s"^). These spectra 
are calibrated using the standard chopper-wheel method which corrects for atmospheric and 
ambient temperature losses to yield the corrected antenna temperature Ta*. 

In the OTF mapping mode, the telescope moves fast and smoothly across the target 
field taking the data continuously, and each map pixel is sampled independently by all 16 
independent detector pixels. Therefore, the use of the OTF observing mode offers a signifi- 
cantly improved calibration, relative pointing accuracy and registration, and a much higher 
dynamic range over the traditional pixel-by-pixel mapping mode. Our OTF observations 
have fully covered the entire stellar disk of each galaxy multiple times. Typically a 10' x 10' 
size box centered on each galaxy is mapped with a scan speed of 45" per second, and the 
data is stored in every 0.25 second. A reference spectrum is obtained after every or ev- 
ery other row of scan at a location 30' away in the azimuth direction. A 6' x 4' region is 
mapped for NGC 4536 because of a telescope problem during the observations. NGC 4567 
and NGC 4568 is an interacting pair observed simultaneously using a single scanning box, 
and one data-cube contains both galaxies. 

The pointing and focus of the telescope are measured at 2-4 hour intervals by observing 
the 86 GHz SiO maser in R-Leo. The measured rms pointing error is ~ 3" in both azimuth 
and elevation, or ~ 5" total. 



3.2. Data Reduction 



The data redu ction is carried out in two steps, initially using the revised -OTFTOOL 

Chung et al.ll2006h and later the GIPS^ (Groningen Image Processing SYstem; Ivan der Hulst et al, 



package. The revised-OTFTOOL reads in the raw OTF data and produces a map 



^http: / / astro. rug. nl/^gipsy 
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Fig. 1. — The sample properties as a function of (a) morphological type (LEDA), (b) major 
axis diameter at 25th mag arcsec"^ in the 5-band (LEDA), (c) corrected total S-band 
magnitude (LEDA), and (d) angular distance from M87. The white region represents the 
28 Virgo galaxies observed, and the shaded region corresponds to the 20 galaxies that are 
detected in CO. 
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after the initial editing and calibration. The resulting data cubes are written out in FITS 
format and are imported to the GIPSY environment for further data reduction and analysis. 

The revised-OTFTOOL is a newly developed program based on the OTFTOOL, which is 
the FCRAO facihty pipehne software for the SEQUOIA OTF data. The OTFTOOL was de- 
signed primarily for the reduction of narrow Galactic emission line data taken with the digital 
backend, and the default data filter is not well suited for the baseline removal in the presence 
of weak, broad emission lines seen in extragalactic CO data. The revised-OTFTOOL includes 
several new functions that arc specifically designed to produce noise-limited output images, 
with improved data filtering and baseline fitting. An improved filtering algorithm in the 
revised-OTFTOOL identifies and removes bad spectra using the rms level, antenna trajec- 
tory, elevation, and system temperature/gain (Tgys), and data containing spikes are identified 
and excluded. A second major improvement is the implementation of a new self-referencing 
method. Rather than using the conventional "OFF" spectrum, our self-referencing method 
constructs the best OFF spectrum from the OTF data itself by choosing the line-free regions 
in the spatial and spectral domain (see Chung ct al. 2005a & 2005b for more detail). This 
self-referencing method produces an OFF spectrum temporally much closer to the ON spec- 
tra, significantly reducing the infiuence of any residual gain changes and thus a significantly 
improved baseline behavior. Finally, the data are normally weighted and mapped onto a 
regular 15" grid. 

The GIPSY package is used to produce the final data cubes from the individual scan 
maps. Occasional bad filter bank channels arc removed and replaced with new data generated 
by interpolating two or more adjacent channels. Interpolation should produce a reasonable 
result since the channel separation {AV ~ 13 km s~^) is small compared with the intrinsic 
CO line width in these galaxies. Data cubes including only the CO emission, used for further 
spectral analysis, are created in two steps. First, a data mask for the CO emission region 
is created through an iterative algorithm that identifies the signal regions with criteria of 
1.5 rms and minimizes the noise in the line- free regions through a removal of a low-order 
(n < 2) polynomial baseline in the frequency domain. Then the final "signal-only" data cube 
is obtained from the noise-minimized data cube by excluding the noise part and retaining 
the emission regions. For 6 galaxies which have large inclination (> 80 degree) or weak emis- 
sion(S/N ratio < 2.5), Position-velocity diagram (PVD) is used to obtain physical quantities. 



Table 2. Information of CO Atlas 
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Group II** - Figure 15-17 
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Note. — ^ Position-velocity diagram (PVD) of Group I galaxy is shown the central slice 
of the dataset along the major axis, and that of Group II galaxy is the integrated PVD. The 
integration is done over some minor axis length for each galaxy (1 or 2 beam- width for highly 
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inclined galaxies and up to 7 beam- width for weak CO and face-on galaxies). 
* NGC 4567 and NGC 4568 are shown in the same map. 

** Group II galaxies have weak CO emission and do not show much emission in the channel 
maps. Therefore only an integrated PVD is shown. 



4. RESULTS 



Among the 28 galaxies observed, 20 galaxies including a galaxy pair of NGC 4567 and 
NGC 4568 are detected in CO emission. The CO detected galaxies are classified into two 
groups : (a) 14 galaxies with strong emission features in the channel map with S/N ratio 
> 2.5 (Group I); and (b) 6 galaxies with large inchnation > 80 degree or weak CO emission 
with S/N ratio < 2.5 (Group II). We have produced a CO atlas for these galaxies and derived 
their CO properties. NGC 4567 and NGC 4568 are an interacting pair, poorly resolved by 
our spatial and velocity resolutions. Therefore, their measured and derived properties are 
reported as a single object. The results of NGC 4536 are derived from 6' x 4' size data cube. 
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Fig. 2. — Group I-l. CO atlas and channel maps of NGC 4254. CO intensity map (top 
left), velocity field map (top right), CO contours overlapped on the optical S-band image 
(middle left), position- velocity diagram which is the central slice along the major axis of 
data cube (middle right), global CO line profile (bottom left), and radial CO distribution 
(bottom right). The line of radial profile is the average of the east and west-side values for 
each radial position. The channel maps are shown over the velocity range of emission. The 
numbers on top right of the first and second channels represent the velocities in km s~^, and 
their difference is the channel velocity separation. 
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Fig. 3. — Group 1-2. Same as Figure 2 for NGC 4302. In the optical image, the companion 
galaxy NGC 4298 is visible in the east. 
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Fig. 4.— Group 1-3. Same as Figure 2 for NGC 4303. 
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Fig. 5.— Group 1-4. Same as Figure 2 for NGC 4321. 
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Fig. 6.— 



Group 1-5. Same as Figure 2 for NGC 4501. 



- 19 - 




(o'ooos) '^^a 




(O'OOOE) 'aaa 



(o'ooos) '^sa 



(iff) ^itsuap XTiij 



Fig. 7.— Group 1-6. Same as Figure 2 for NGC 4527. 
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Fig. 8.— 



Group 1-7. Same as Figure 2 for NGC 4535. 



- 21 - 



o 

o o 



—1 W 




- . - 

■V, 
,1* -■ ■ 

"i 1 1 1 1 1 1 1 1 1 1 1 1 r 


- - 

" 'ts. 

r 1 1 1 1 1 1 1 1 r 


fi 1 1 1 1 1 1 1 1 1 1 1 1 r 


■ & . 

fi 1 1 1 1 1 1 1 1 1 1 1 1 r 


^ " ■ 

■ ' < 

fi 1 1 1 1 1 1 1 1 1 1 1 1 r 


P ■ _ 

"i 1 1 1 1 1 1 1 1 1 1 1 1 r 


■ ? 

^ ■ 

[i 1 1 1 1 1 1 1 1 1 1 1 1 r 


[i 1 1 1 1 1 1 1 1 1 1 1 1 r 


fi 1 1 1 1 1 1 1 1 1 1 1 1 r 


fi 1 1 1 1 1 1 1 1 1 1 1 1 r 


- CD 

in o 

"i 1 1 1 1 1 1 1 1 1 1 1 1 r 


- - 

fi 1 1 1 1 1 1 1 1 1 1 1 1 r 


^"'.^ 

_ F 

[i 1 1 1 1 1 1 1 1 1 1 1 1 r 


°-<-> 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 


■ » c 

fi 1 1 1 1 1 1 1 1 1 1 1 1 r 


" O ■ A 




- . ^ - 


; S ; 


■ 









1 3 






Li 




er (arcn 










a 














(o'oooe) '^^a 






i # 

I 




□ 







(0 0003) (00003) "^SQ (^f) ^;rsKap XTiy 



Fig. 9.— 



Group 1-8. Same as Figure 2 for NGC 4536. 
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Fig. 10.— Group 1-9. Same as Figure 2 for NGC 4567 & NGC 4568. NGC 4567 is in the 
center, and NGC 4568 is in the south. The PVDs and radial CO distributions are obtained 
according to each galaxy's center and position angle. 
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Fig. 11.— Group I-IO. Same as Figure 2 for NGC 4569. 
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Fig. 12. — Group I-ll. Same as Figure 2 for NGC 4647. In the optical image, the elhptical 
companion galaxy NGC 4649 is visible in the east. 
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Fig. 13.— Group 1-12. Same as Figure 2 for NGC 4654. 
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Fig. 14.— 



Group 1-13. Same as Figure 2 for NGC 4689. 
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Fig. 15.— Group II-l. CO atlas of NGC 4298(left) and NGC 4402(right). Optical 5-band 
image (top left), integrated position- velocity diagram (top right), global CO line profile 
(bottom left), and radial CO distribution (bottom right). In the optical image of NGC 4298, 
its companion galaxy NGC 4302 is visible in the west. 
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Fig. 16.— Group II-2. Same as Figure 15 for NGC 4419(left) and NGC 4438(right). 
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4.1. CO Atlas: Map Descriptions 



Figures 2-17 are the CO atlas and selected channel maps of Group I and Group II 
galaxies obtained from our OTF mapping observations. The atlas figure for each Group I 
galaxy consists of 6 maps: a CO intensity map (top left), a velocity field map (top right), 
an optical image with CO contours (middle left), a position- velocity diagram derived along 
the major axis of the stellar disk (middle right), a global CO line profile (bottom left), and a 
radial CO profile (bottom right). Selected channel maps are also shown for Group I galaxies. 
For Group II galaxies, CO emission is too faint to yield an intensity map and velocity-field 
map. Therefore, only an optical 5-band image from NASA/IPAC Extragalactic Database 



(NEDB; top left), a position- velocity diagram derived by integrating along the minor axis 
(top right), a global CO line profile (bottom left), and a radial CO profile (bottom right) 
are shown. 

Each map shown (including the channel maps) has a dimension of 6' x 6' (4' x 4' for 
NGC 4536) since no emission was detected outside of this region for all galaxies. The veloc- 
ity range shown for the global CO profile and channel map is about 500 km s~^, except for 
NGC 4501 which has a CO line width of about 520 km s~^. The contour levels shown are 
listed in Table 2. 

CO intensity map: a CO intensity map is produced by summing all emission features 
in the channel maps. The contour levels shown are listed in K km s~^, and the CO line 
intensity is in the T^ scale. 

Velocity field map: a mean velocity field map is derived through single Gaussian fit to 
the observed line profiles. 

Optical image: An optical 5-band image is obtained from NED. 

Position-velocity diagram: a PVD is extracted from a row of pixels along the major axis, 
i.e., a central slice of 15 arcsec width is shown in ATLAS, for Group I galaxies. For Group 
II galaxies, whose inclination is too high or whose CO emission is too weak in the individ- 
ual channel maps, an integrated PVD is produced by summing the data along the minor 
axis. The integration is done over some minor axis length for each galaxy (1 or 2 beam- 



•^http: / /nedwww. ipac.caltech.edu 
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Fig. 17.— 



Group II-3. Same as Figure 15 for NGC 4548(left) and NGC 4579(right). 
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width for highly inchned galaxies and up to 7 beam- width for weak CO and face-on galaxies). 



Global CO line profile: For Group 1 galaxies, a global CO line profile is extracted from 
the data cube. For the Group 11 galaxies, CO flux density is integrated manually along 
the minor axis of the PVD , and the global line profile is computed (e.g., iBroeils fc Rhee 
19971 : iRhee &: Broeilsll2005l ). The uncertainty in the integrated line flux in each channel is 
calculated as 



O": 



,cal 



0"; 



T* X 0.1, 

Trms 
i -I 



„bsln 



(T* - Tio)2 + (T* - Ta)2 + (T* - T,^) 



(2) 
(3) 

(4) 



where 0"*°* is the total uncertainty, cr?**', o"[™^, and af^^'^ are the calibration uncertainty, 
thermal noise in the channel maps, and the baseline subtraction noise. F or the calibration 
uncert ainty, we assign a 10% absolute uncertainty for all spectra following iKenney fc Young 
( 1l988l ) (Eq. [2]). T* is CO line intensity, T[™** is map noise level, and Tjo, Tji, and Ti2 are 
CO intensities o f the i-th channel deriv ed by the baseline subtraction of 0, 1, and 2 order, 
respectively (see lKenney fc Young] Il988l ). 



Radial CO profile: deriving a radial surface density proflle using a concentr i c ellip se fltting 
method does not work well for edge -on or highly inclined galaxies. IWarmeld (119881 ) applied 
Lucy iteration method (ILucyl Il974l ) to obtain HI ra dial surface den sity profile of edge on 
galaxies, and we derived radial CO profiles following IWarmelsl (119881 ). 



Channel map: The velocity range of each channel map is about 500 km s ^ (740 km s 
for NGC 4501), which is the same as that of the global CO line profile. 



4.2. CO Properties 

CO properties for the 20 CO detected Virgo spiral galaxies are summarized in Table 3. 
The column entries are: (1) NGC number; (2) rms noise measured from emission-free re- 
gions; (3) line width measured at 20% level; (4) line width at 50% level; (5) systemic velocity 
derived from the CO velocity profile; (6) total CO line fiux; (7) molecular hydrogen mass; 
(8) effective CO diameter; and (9) isophotal CO diameter. 
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Total CO flux & H2 mass 



To convert the CO intensity in T^ into the flux density unit, a cahb ration factor of 42 Jy K^^ 
obtained using the FCRAO 14-m telescope (IKenney fc Young|l988l ) is appl ied. We derive H2 
mass es from total CO flux measured assuming a linear conversion relation (IKenney fc Young 
1989h . 



H2 



3.9 X lO-iV'Sco(M0) 



(5) 



where x is the conversion factor and d is dist ance to the source in ni egaparsec (Mpc). We 



adopt X = 3 X 10^0 cm'^ ( K[Tr] km s"^ flYoung fc Scovillelll99lh. The distance to the 



Virgo cluster is somewhat uncertain because of the large d epth effect (|Y"asuda et al.lll997h and 



is generally thought to b 


1999 




Sofue et al. 


2003: 



Young et al. 



1995 : 



Sakamoto et al. 



of 20 Mpc for an easier comparison with the FCRAO Extragalactic CO Survey (see § 14.31) . 



CO line width 

CO line width is measure d at 20% and 50% level of the line peak, on each side of the line, 
following the deflnition of iRhed (119961 ). The line proflle is divided into two equal velocity 
bins, and the peak fluxes "T^^^ and Tj^fg^ are determined separately on each half of the line 
proflle. The 20% and 50% velocities {y"^^'^" and V^'^'^°) represent the velocities at which the 
line proflle reaches the 20% and 50% of the peak value on the respective high and low velocity 
side, approaching from the line edge to the center of the line proflle. A linear interpolation 
procedure is used in this calculation. 



The flnal value of line width of each cutoff-level is determined as 



W 



obs 
20 



Y20% 
^ high 



V 



20% 
low 



Y50% _ y50% 
high low ■ 



(6) 
(7) 



The uncertainty in the line width is est imated with 1 a uncertain ty in each case where 
cr is a mean rms noise of line-free channels (iRhee &: van Albadalll996l ). 



The line width is corrected for instrumental broadening following the method described 
by IVerheijenI (119971 ). The correction for instrumental broadening (in km s~^) is computed as 
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W20 = W°Jj^ - 35. 



W50 = W^^' - 23.5 




(9) 



where AV^ is the velocity resolution in km s ^ . 



CO systemic velocity 



Again, following the procedure described by IVerheijenI (119971 ). we compute the systemic 
velocity of each galaxy as 



V, 



sys 



l^low "T Vhigh+ ^low + V high J 



(10) 



CO diameters 



To describe CO extents, we derived effective diameters and isophotal diameters. Effective 
CO diameter, D^ , is defined as a diameter which encloses 70% of the total emission (see 



Young et al.lll995l ). The CO isophotal diameter, DJ^g, is defined as a diameter where the 
mean face-on surface density of H2 falls to 1 MqPc"^ which corresponds to / Tr dv = 0.21 

? cm 

(K[Tr] km s-i )-i dYoung fc Scovillelll99lh . 



K km s ^, or 6.3x 10^^ cm ^ for a CO -H^ proportionality factor of x 



3x 10^0 cm-2 
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4.3. Comparison With Other Observations 



Previous single dish and interferometric CO measurements exist for many of our target 
galaxies. A compa rison of our results with those of the FCRA O Extragalactic CO Survey 
dVoung et all llQQsh and the BIMA SONG jHelfer et all WM) provides an important test 
of potentially important systematics associated with these two widely utilized extragalactic 
CO surveys. 



4-.3.1. Comparison With the FCRAO Extragalactic CO Survey 



The FCRAO extragalactic CO survey (lYoung et al.lll995l ) data on 300 external galaxies 



is the most extensive and most widely referenced database for CO emission measurements 
in 300 external galaxies. A comparison with our CO imaging survey of Virgo cluster spi- 
rals provides an important verification of this earlier survey, which employed sampling and 
modeling rather than full imaging, and offers a quantitative constraint on the systematics 
resulting from the emission modeling used by Young et al. in estimating the total CO lu- 
minosity. Since both surveys used the same telescope and calibration method as well as 
adopting the identical distance to the Virgo Cluster, any hardware-dependent systematics 
are minimized. Therefore any measured difference is narrowly constrained to the differences 
in the observing methods and the modeling used by Young et al. NGC 4536 is excluded in 
this comparison because the small area coverage of our maps may have adversely affected 
the baseline subtraction process. 



Table 3. Derived CO Properties 



Galaxy rms Wm^ V.y,'' Sco A/h^ Df^'' Dg^" 

(K) (km s^-^) (km s^^) (km s"-"-) (Jy km s^-"-) (10**Mq) (arcmin) (arcmin) 



Group I 



NGC 4254 


0.016 


221±1 


209±4 


2392 


2830±290 


134±14 


2.6 


5.1 


NGC 4302 


0.013 


353±1 


344±4 


1147 


540±150 


25±7 


2.7 


3.6 


NGC 4303 


0.018 


162±2 


145±4 


1566 


1920±200 


90±9 


2.2 


3.9 


NGC 4321 


0.021 


239±2 


217±4 


1580 


2390±270 


113±13 


1.8 


4.5 


NGC 4501 


0.017 


518±2 


502±3 


2276 


2130±240 


lOlitll 


2.5 


4.8 


NGC 4527 


0.017 


376±2 


358±5 


1736 


1260±160 


60±7 


1.7 


3.4 


NGC 4535 


0.029 


252±5 


235±5 


1960 


1450±250 


68±12 


2.3 


3.9 


NGC 4536* 


0.010 


318±2 


267±5 


1802 


390±50 


19±3 


0.8 


1.9 


NGC 4567" 


0.018 


297±3 


272±4 


2242 


1500±170 


71±8 






NGC 4569 


0.022 


311±1 


297±5 


-219 


900±130 


43±6 


1.3 


2.9 


NGC 4647 


0.016 


156±1 


145±12 


1426 


620±100 


29±5 


1.2 


2.8 


NGC 4654 


0.014 


285±1 


270±5 


1039 


710±90 


33±4 


1.8 


3.4 


NGC 4689 


0.020 


182±1 


171±3 


1621 


510±80 


24±4 


1.6 


2.8 


Group II 


















NGC 4298 


0.012 


270±4 


214±10 


1145 


410±60 


20±3 


2.1 


3.6 


NGC 4402 


0.025 


267±3 


251±5 


242 


520±100 


24±5 


1.4 


1.5 


NGC 4419 


0.016 


318±2 


290±13 


-218 


390±100 


18±5 


0.9 


2.0 


NGC 4438 


0.018 


272±5 


151±8 


161 


300±70 


14±3 


1.1 


2.3 


NGC 4548 


0.015 


234±6 


217±5 


484 


710±90 


34±4 


2.1 


3.7 


NGC 4579 


0.014 


329±1 


317±5 


1502 


500±70 


23±3 


2.1 


3.1 



Note. — Linewidths are corrected for the instrumental broadening (see Eqs. |8] &; [SJ. Uncertainties are 
derived by following Rhee & van Albada (1996). 
CO systemic velocity 

Diameter which contains 70% of the total CO flux 
" Diameter where the mean face-on surface density of H2 falls to 1 Mq pc"'^. H2 surface density distribution (in 
unit of M0 pc~^) is derived from CO radial distribution (in unit of K km s~^) assuming a constant CO to H2 
conversion factor within a whole galaxy (x = 3 X 10^° ^^-2 ( K[Ta] km s"! )-i l lYoune fc ScovillellTggih ) . 
* CO properties of NGC 4536 is derived from the datacube mapped by 6' X 4' size. 

**The entries for NGC 4567 are the total values for NGC 4567 and NGC 4568 as a pair. CO diameters are not 
given for these galaxies. 



Table 4. Comparison with 



Young et al, 



(119951) 



Galaxy FCRAO 14-m (OTF) " FCRAO 14-m (Position Switching) 



cjOTP nOTF nOTP oPS.flt „„r oPS.obs nPS-Ht nPS-flt 

(Jy km s^-"-) (arcmin) (arcmin) (Jykms^-') (%) (Jy km s"-"-) (arcmin) (arcmin) (%) 





(1) 


f3"l 




('5'l 


y") 


K ' ) 








Group I 




















NGC 4254 


2830±290 


5.4 


2.6 


3000±850 


37 


1110 


5.3 


3.4 


6 


NGC 4302 


540±150 


4.6 


2.7 


620±100 


100 


620 


2.0 


1.8 


15 


NGC 4303 


1920±200 


4.2 


2.2 


2280±470 


50 


1140 


3.4 


2.2 


19 


NGC 4321 


2390±270 


5.4 


1.8 


3340±920 


39 


1300 


4.5 


3.4 


40 


NGC 4501 


2130±240 


5.2 


2.5 


2220±480 


59 


1310 


3.8 


2.9 


4 


NGC 4527 


1260±160 


3.8 


1.7 


1800±410 


92 


1660 


3.4 


2.0 


43 


NGC 4535 


1450±250 


4.1 


2.3 


1570±410 


38 


600 


3.9 


3.5 


8 


NGC 4567 * 


1500±170 






1550±210 


81 


1250 






3 


NGC 4569 


900±130 


3.1 


1.3 


1500±260 


82 


1230 


3.2 


1.8 


67 


NGC 4647 


620±100 


3.2 


1.2 


600±120 


66 


400 


2.2 


1.6 


-3 


NGC 4654 


710±90 


3.6 


1.8 


730±150 


62 


450 


2.6 


2.4 


3 


NGC 4689 


510±80 


3.0 


1.6 


710±150 


50 


360 


2.4 


2.4 


39 


Group II 




















NGC 4298 


41Q±60 


1.5 


2.1 


660±110 


86 


570 


2.4 


1.8 


61 


NGC 4402 


520±100 


1.8 


1.4 


630±120 


98 


620 


2.4 


1.4 


21 


NGC 4419 


390±100 


1.5 


0.9 


920±190 


91 


840 


2.2 


1.4 


136 


NGC 4438 


300±70 


1.4 


1.1 


210±40 


94 


200 


1.3 


1.2 


-27 


NGC 4548 


710±90 


2.4 


2.1 


540±140 


60 


320 


3.3 


3.0 


-24 


NGC 4579 


500±70 


2.5 


2.1 


910±200 


47 


430 


2.7 


2.6 


82 



00 



Note. — " Our new OTF observation results 
' Position-Switching observation results by Young ct al. (1995) 
Col. (1) — NGC number 

Col. (2) — Total CO flux, S^q^, and its rms, ctotFi derived from our OTF observations 

Col. (3) & (4) — The isophotal and effective CO diameter, DPJI" and 0°/^, derived from our OTF observation. To compare with 
Young et al. (1995), Dj^^^ in here is derived at the diameter where the face-on CO integrated intensity falls to 1 K(T^) km s~^. 
Col. (5) — The fitted total CO flux, Sg^^', and its rms, a^t, by Young et al. (1995) 



Col. (6) — Scale factor which is the percentage of the total CO emission sampled by the observations by Young et al. (1995) 
Col. (7) — Recomputed Young ct al. PS observed CO flux, S^l,"''^, using Col. (5) and Col. (6). 

Col. (8) &: (9) — The isophotal and effective CO diameter, D?^ *"' and D^J '^*, which is the diameter where the face— on CO integrated 
intensity falls to 1 K(T^) km and the diameter which contains 70% of the total CO flux for the best-fitting model, respectively 
(Young et al. 1995). 

Col. (10) — Fractional ratio of the difference between S^o^ and Sg^,-***, calculated as (Sgf,''* - Sg5^)/Sg5^ X 100. 

* NGC 4567 and NGC 4568 cannot be separated in our OTF observations, and the summed total flux is reported. For the PS results, 

the equivalent summed CO flux of NGC 4567 and NGC 4568 is listed (Young et al. 1995). 
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A comparison of the FCRAO 14-m OTF and PS measurements are summarized in Ta- 
ble 4. Both the raw PS (S^o°''^) and th e model fit (Snn^^) results are listed, using the scale 
factor (scf) given in the original paper (lYoung et al.lll995l ). An isophotal diameter D^J^ is 
defined as the diameter where the face-on CO integrated intensity falls to 1 K(T^) km 
while an effective diameter is defined as the diameter which contains 70% of the total 

CO fiux, following the definition by Young et al.(1995). The fractional ratio between the 
OTF and PS measurements, r, is calculated and presented to show the difference between 
the two methods of the total CO flux. Comparisons of the OTF and PS measurements are 
also shown in Figure | 



The panel (A) in Figure [TS] shows that the observed PS line fluxes are systematically 
smaller than the OTF measurements, but this is expected since "observed PS" data include 
only partial measurements taken along the major axis. A better agreement is expected for 
galaxies with a high inclination and a small size. And indeed galaxies with a smaller D25 
and a larger inclination appears to have a smaller discrepancy between S^o°^'' and S^g^. 
The model-flt measurements shown in panel (B) are more consistent with the OTF measure- 
ments. The OTF and PS model-flt measurements of 9 galaxies are consistent with the total 
CO flux to within la. Out of the remaining 9 galaxies, 7 show larger Sqq'^* than Sqq^ by a 
factor of 1.4 - 2.4. And, the last two galaxies (NGC 4438 and NGC 4548) have larger SgJ'' 
than Spo'^* by a factor of 1.4. Among the 9 galaxies for which the discrepancies are larger 
than la, there are 4 galaxies which disagree by > 2a uncertainty, and the 4 galaxies have 
larger Sqq^^ than S^q^ by a factor of 1.6 - 2.4. We predict that peculiar CO distributions 
such as ring- or bar-like structures (e.g.. Young et al. 1995) can affect the model fltting 
process, and more analysis will be done in another paper (Chung et al. 2010). 
Panels (C) and (D) in Figure [18] show the comparisons of CO diameters derived from the 
OTF and PS observations. For Group I galaxies detected with high S/N ratio, isophotal di- 
ameters derived from the OTF data (D^J^) are larger than those of the PS data (D?^-^*). In 
contrast, the effective diameters of OTF (D^^) are smaller than DgJ '^*. An explanation for 
these apparently puzzling trends is that CO emission is spread over a larger extent but has a 
more concen t rated central component in the OTF maps when compared with the results of 
Young et al.l (119951 ). The effective diameters are smaller than 3 arcmin in every case. 



and the isophotal diameters are located between 1/2 and 1 of the optical diameter D25 in 
most cases. D^J^ for Group II galaxies are systematically smaller than D?^'^*, and this may 
be the result of low S/N ratios of the OTF data. 
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Fig. 18. — Comparison of the data with the FCRAO extragalactic CO survey (jYoung et al. 



19951 ). Here, "PS" denotes the position-switch measurements reported in the FCRAO extra- 
galactic CO survey data while "OTF" refers to our On-the-Fly data. Panel (A) compares the 
total CO fluxes of the OTF measurements with the PS measurements. Panel (B) compares 
the total CO fluxes of the OTF measurements with the model fit corrected PS data. Panels 
(C) and (D) compare the isophotal and effective diameters derived from the OTF and the PS 
data, respectively. Isophotal diameters D^I^ D^^-^* are defined as the diameter where 
the face-on CO integrated intensity falls to 1 K(T^) km s~^. Effective diameters and 



Dgg'^* are the diameter which contains 70% of the total CO flux. The dotted line represents 
a perfect linear relationship. Filled and open symbols denote Group I and Group II galaxies, 
respectively. 
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4-3.2. Comparison With the BIMA SONG Survey 



Three galaxies imaged in the OTF mod e (NGC 4303 NGC 4321, & NGC 4569) were 
also observed by the BIMA SONG project Jflelfer et aljl2003h . The total BIMA SONG 
CO fluxes reported for NGC 4303, NGC 4321, and NGC 4569 are 2427 ± 145, 2972 ± 319, 
and 1096 ± 137 Jy km s~^, respectively. These integrated CO line fluxes are systematically 
slightly larger (~ 25%) compared with our OTF measurements and are between the OTF 
and PS measurements on average. The apparent difference may be rooted in the systematic 
difference in calibration. Because of the small number of objects in common, we can only 
report a broad agreement. 



5. DISCUSSION 

5.1. Notes on Individual Galaxies 

In this section, we offer supplementary information on the individual galaxies and offer 
comparisons with other CO, HI, and Ha observations. In particular, we are interested in 
the environmental influence on the molecular ISM, and we make notes of any peculiar CO 
distribution that may be related to the cluster environment. 



5.1.1. Group I Galaxies 



NGC 4254: This Sc galaxy is located 3.3° northwest of M87. ICavatte et al.l (119901 ) re- 
ported that it has a very asymmetric HI distribution with a sharp edge on the side pointing 
toward M87 and an extended HI feature pointing away from M87, as if it were compressed 
by the interga l actic medium. In our CO intensity map, CO appears to extend to the south. 
Chemin et al.l (120061 ) showed an asymmetric structure in Ha velocity field. Streaming mo- 
tions along the spiral arms and steep velocity rise are observed in their atlas. Our CO 
velocity field map and position-velocity diagram are in good agreement with the Ha kine- 
matics. Interferometric high resolution imaging of the central region shows clumped CO 
distribution along the spiral arms, a bar - like elongation fea ture, and slight depression at the 
dynamical center (ISakamoto et al.lll999l : ISofue et al.ll2003l ). 



NGC 4302: This Sc galaxy is a dusty edge-on spiral located 3.1° northwest of M87. Its 
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CO distrib ution in our CO int ensity map consists of two molecular peaks with a central 
depression. lYoung et al.l (Il995l ) modeled the CO distribution as being offset by 0'.45 from 
its center. Our new data suggest that this earlier study might have detected only one of the 
two peaks. 



NGC 4303: This Sc spiral is located 8.2° south of M87. ICavatte et al.l fll990[ ) reported a 
symmetric HI distribution with a central depression. Our atlas shows that CO emission is 
distributed along the spiral arms, an d a bar-like structure is displaced in position angle by 
about 20 degrees. iHelfer et al.l (|2003l ) reached the same conclusion using the BIMA SONG 
survey data. The CO emission is concentrated in the central region, unlike the HI (see Cay- 
atte et al. 1990), but the HI and CO velocity fields are in good agreement. A recent burst of 
st ar formation caused b y tidal interactions with two nearby companions has been suggested 
bv lCavatte et al.l Jl99oh . 



NGC 4321: This Sc galaxy is located 3.9° north of M87. Its HI emission is distributed 
along the spiral arms, and the oval dist ortion in the inner p art due to the presence of a 
bar is shown from the kinematic pattern (jCayatte et al.lll990l ). Our CO image shows a well 
confined CO distribution on the stellar disk, plus an asymmetric extension to the southern 
spiral arm. It also has an elongated bar like structure along the major axis. The CO veloc - 
ity field shows an oval distortion in the inner part, similar to the HI (jCayatte et al.l 119901). 



The BIMA SONG survey detected CO at the center and along the sp iral arms (jHelfer et al. 



20031 ) . in a good agreement with our results. ISakamoto et al.l (119991 ) reporte d a prominent 



pair o f nuclear CO arms and a sharp condensation of CO at the nucleus. IChemin et al. 



(120061 ) found three different pattern speeds, which is a sign of a streaming motion. 



NGC 4501: This Sbc spiral is located 2.1° north of M87. ICavatte et all (Il990f ) reported 
that this galaxy points its steep HI edge toward M87, similar to NGC 4254, and suggested 
these two galaxies as examples of enhanced star formation c aused by compression o f the ISM. 
Our CO intensity map shows a well -centered dis t ribut ion. ISakamoto et al.l (|l999l ) reported 
a concentrated CO morphology, and lSofue et al.l (|2003l ) found that the north-eastern arm is 
much brighter in CO than the south-western arm. The CO velocity width of NGC 4501 is 
one of the largest among the galaxies observed. The Ha v elocity field appears regular, but 



the position- velocity diagram shows a complex kinematics (IChemin et al. 



20061 ). 



NGC 4527: This Sb galaxy is located 9.8° south of M87. iKennev fc Yound (119881 ) reported 
that it has a uniform disk component. Our OTF map is consistent with the Young et al. 
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data, but the OTF map shows that there is a central CO concentration along with a bar 
along the major axis and an asymmetric extension to the southwest. 

NGC 4535: This SBc spiral is locat ed 4.3° south of M87 . This galaxy shows an undisturbed 
HI distribution with a central hole (jCayatte et al.lll990l ). iHelfer et al.l (12003 ') reported that 
CO is distributed at the center and along the spiral arms. Isofue et ah (2003) classified this 
galaxy as a typical single-peak type with offset bars. Our OTF data is noisy, and the signal 
is very weak. The CO intensity map shows a bar-like elongated structure in the central re- 
gion, and the p ositio n- velocity diagram shows evidence for highly disturbed gas kinematics. 
Chemin et al.l (120061 ) reported a perturbed velocity field and streaming motions along the 



arms. 



NGC 4536: This Sc galaxy is located 10.2° south of M87. Our map covers only a 6' x 4' 
region, but the CO emis sion appears to be fully covered when compared with the CO map of 
Kenney fc Young! ( 1l988l ). CO emis sion is strongly co ncentrated on the galactic center, and 
this agrees well with the results by lSofue et al.l (120031 ). 



NGC 4567 & 4568: Thes e two Sc galaxie s are l ocated 1.8° , and they are not well separated 
spatially or kinematically. ICayatte et al.l (jl990l ) suggested that the HI emission displaced 
tow ard the south of NGC 4 567 c ould be a sign of a tidal interaction between the two galax- 
ies. iKoopmann fc Kenneyl (120041 ) also suggested ram pressure effects and a tidal interaction 



between the two. However , Ha velocity field does not show any clear signs of velocity 
disturbanc e s (jChemin et al.ll2006l ). Higher angular resolution HI and CO observations by 
lono et al.l (120051 ) also show that the inner gas disks show little signs of tidal disturbance, 
with a symmetric bar-like or spiral-like features. Unlike many other observed Virgo galax- 
ies, they both show CO emission extending out to the outer optical radii in our OTF map, 
including where the two disks overlap. 



NGC 4569: This Sab spiral is an Hl -anemic galaxy (Ivan den Berghl Il976l ) located 1.7° 
northeast of M87. ICayatte et al.l (Il990l ) reported that its HI disk is severely stripped. In- 
terferometric C O imaging found a symmetric CO di stribution with two peaks and a cen- 
tral depression (iHelfer et al.l l2003l : ISofue et al.l l2003l ). Higher resolution CO imaging by 
Nakanishi et al.l (120051 ) found a highly concentrated CO distribution in the circum-nuclear 
region with two off-center peaks. Our CO atlas shows a centrally concentrated distribution 
with a much smaller CO radius than the optical radius and a bar-like elongated structure 
across the optical major axis. Its Ha velocity field is perturbed and shows evidence for 
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streaming motions (jChemin et al.ll2006l ). Our CO velocity field map shows that central ve- 
locity is shifted toward the south. 



NGC 4647: This Sc galaxy is located 3.2° east of M87 and has a companion elliptical 
galaxy NGC 4649. N GC 4647 shows a slightly extended HI distribution toward NGC 4649 
(ICayatte et al.lll990l ). However, CO appears to point away from its elliptical companion 
NGC 4649 in our intensity map. The CO position-velocity diagram also appears disturbed. 
The companion galaxy NGC 4649 is undetected in CO. 



NGC 4654: This SBc galaxy is located 3.3° northeast of M87. The HI image shows a 
sharp cutoff on the northwest side and an eastward extens io n, and enhanced s t ar for mation 
activity is also seen in the northwest (ICayatte et al.lll990l ). iKenney &: Young (Il988l ) found 
an asymmetric CO distribution - the position 45" northwest of the nucleus shows a stronger 
line than the central 45", and t his region also displays a peak in the HI, radio continuum, 
and the Ha. ISofue et al.l (120031 ) found a lopsided CO distribution in the inner region and 
suggested that even the nuclear region suffers from the ram pressure effects. Our OTF CO 
maps also shows a clear asymmetry and lopsidedness, but the CO peak is well centered on 
the stellar disk. The lopsided CO extension is both to the northwest and to the southwe st, 
in the direction of M87. The Ha velocity field is strongly perturbed (jChemin et al.ll2006l ). 



NGC 4 689: This Sc galax y is located 4.3° northeast of M87. It has a large, ex t ended 
HI disk (ICayatte et al.lll990l ). High angular resolution CO imaging by ISofue et al.l (120031 ) 
found a lopsided, amorphous CO morphology and without a central peak. Our CO map 



shows a CO peak offset from the optical center and extended to the south. IChemin et al. 



( 120061 ) reported that Ha velocity field is slightly perturbed due to streaming associated with 
pseudo-spiral arms. 



5.1.2. Group II Galaxies 

NGC 4298: This Sc galaxy is located 3.2° northwest of M87 and has a companion NGC 4302 
only 2.3 arcmin away. Stellar asymmetry due to a recent tidal intera ction is seen in optical 
imag es , and a truncat ed gas disk due to ram-pressure is also suggested (IKoopmann &: Kenney 
20041 ). IChemin et al.l (j2006) suggested that its mildly perturbed velocity field may indicate 
a streaming motion or a locally warped arm. Our PVD shows stronger emission on the east 
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(receding) side than on the west (approaching) side, and CO emission is probably extended 
toward southeast side, in the direction of its companion galaxy NGC 4302. This feature is 
in good agreement with the Ha kinematics discussed above. 



NGC 4402: This edge-on Sc galaxy is loca ted 1.4° northwest o f M87. Although the high 
velocity part of the HI emission is missed by lCayatte et al.l (1l990l ). it is obvious that the HI 
disk of this galaxy is significantly smaller in extent than that of the optical d isk. Our CO 



positi on- velocity diagram shows that CO emission is asymmetric to one side, and I Young et al. 



(I1995I ) also suggested the distribution model to be offset by 0'.20 from the center. A 10" 
radius nuclear disk and a more extend ed molecular disk with a 30" diameter are found by 
higher angular resolution observations (ISofue et al.ll2003l ). 



NGC 4419: This Sa galaxy is located 2.8° north of M87. iKennev k Yound (Il988f ) pro- 
posed that NGC 4419 is on a radial orbit which passes very close to the cluster core and 
strongly interacts with the intracluster medium, resulti ng in a large CO/H I flux ratio and 
a significant CO asym metry. CO distrib ution model of lYoung et al.l (119951 ) shows an offset 
0'.25 from center, and ISofue et al.l (120031 ) reported that the outer molecular disk is lopsided 
toward the northwest. Our CO position-velocity diagram shows a compact distribution and 
a rapidly rising rotation speed in the central region. 



NGC 4438: This Sb galaxy is located 1.0° northwest of M87. ICavatte et al.l (119901 1 re- 
ported a highly asymmetric HI distribution with an extension pointing , and a much smaller 
HI disk than the optical dis k. Its ionized gas has an off-plane filamentary morphology to th e 



east and south of the disk (IChemin et al.l l2006l : iKenney et al.l Il995l : iKenney fc Yald |2002|) . 



Our CO position-velocity diagram shows that its molecular center is offset from the optical 
center by ~ 0.5 arcmin to the east. It also shows a highly disturbed structure. 



NGC 4548: This SBb galaxy is located 2.4° northe ast of M87. This galaxy is a nearly face- 
on ba rr ed (SBb) spiral which is severely Hl-deficient (Ivan den Berghlll976l : iGiovanelli &: Haynes 
19831 ). ICayatte et al.l (Il990l ) reported a ring st ructure in the HI distri bution. The BIMA 
SONG survey detected CO only at the center (IHelfer et al.l 120031 ) . and ISofue et al.l (120031 ) 
also reported a very weak and highly concentrated CO distribution. Our CO data is weak 
and too noisy to determine the molecular distribution well. However, its global CO line 
profile shows moderately symmetric double horns and CO appears to have comparatively 
large extent from the integrated PVD and radial profile. 
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NGC 4579: This Sa b galaxy is located 1.8° east of M87 and is another Hl-anemic spiral 
( Ivan den Berghl Il976l ). Its HI distribution shows a ring-like structure, and it is the most 
severely stripped ga laxy with signs of unusual nuclear activity, possibly fueled by gas inflow 



dCavatte et aDll99oh. The BI MA SONG survey detected CO only at the center jHelfer et al. 



20031 ). and lSofue et al. ( 2003 ) reported that elongated CO distribution is displaced from the 



optical bar axis by about 30 degrees. Its Ha velocity field shows evidence for perturbed, 
streaming motions (jChemin et al.ll2006l ). Our CO position- velocity diagram shows a severely 
disturbed morphology and kinematics. 



5.2. Summary of CO distribution 

Here, we summarize the CO distribution of the Virgo spiral galaxies shown by our 
survey. 



1. CO is confined to the galactic center and disk in most galaxies. However, slight asym- 
metric distribution and significant structures as well as kinematic disturbances are 
frequently shown. This is consistent with previous results of CO studies (e.g.. Young 
et al. 1995; Heifer et al. 2003). 

2. Among the 14 galaxies which have CO intensity maps, 11 galaxies have single CO peak 
center (NGC 4254, NGC 4303, NGC 4321, NGC 4501, NGC 4527, NGC 4535, NGC 
4536, NGC 4567, NGC 4568, NGC 4569, and NGC 4654). NGC 4302 has twin peaks 
with a central depression, and NGC 4647 and NGC 4689 have no notable CO peak at 
the center. 

3. Five galaxies appear to have bar-like elongated structure (NGC 4303, NGC 4321, NGC 
4527, NGC 4535, and NGC 4569). Extended CO distributions to one side are also seen 
in 6 galaxies (NGC 4254, NGC 4298, NGC 4527, NGC 4654, NGC 4689, and NGC 
4402). 

4. CO emission in NGC 4298 shows an asymmetric extension toward the southeast, which 
is the direction of its companion galaxy NGC 4302. NGC 4647 shows a CO distribution 
pointing away from its elliptical companion NGC 4649. 

Figure [19] shows the total CO extent and morphology of each of the 14 Group I galaxies 
at their proper position in the Virgo cluster. There are several noteworthy trends: 
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N4536 



13^ 4™ 48=^ 32"^ 12^16"^ 12^ 0™ 

R.A. (2000.0) 

Fig. 19. — Spatial distribution of the 14 Group I galaxies in the Virgo Cluster. Individual 
galaxies are 20 times enlarged to be shown more clearly. Contour levels correspond to a H2 
column density of 0.75, 1.50, 2.25, and 3.00 10^° cm~^. The radius of each dashed circle is 
3, 5, and 10 degrees, respectively. 
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1. It appears that CO molecules are not strongly affected by the cluster environments 
such as distance from the cluster center of M87. 



2. Galaxies located in the western side of Vi rgo Cluster are kno wn to have a larger HI 
disk than the galaxies in the eastern side (jCayatte et alJIlOQCl ). and we find the same 
trend in the CO disk size as well. 

3. NGC 4298 and NGC 4647 have companion galaxies, but CO emission is extended 
toward and on the opposite to its companion, respectively. Figure [19] suggests that 
both of their extensions point toward the cluster center near M87. NGC 4654 also 
shows a CO extension toward M87. 



6. SUMMARY & CONCLUDING REMARKS 



We have carried out a ^^CO(J=1-0) OTF mapping survey of 28 Virgo cluster spiral 
galaxies. Although the importance of molecular contents in galaxy evolution is widely rec- 
ognized, systematic CO imaging surveys of a large galaxy sample covering the full extent of 
the stellar disks are still rare. Taking advantages of the OTF mapping mode of the FCRAO 
14-m telescope, we have imaged a large area of each galaxy covering the entire stellar disks 
in a relatively short period of time compared with the more traditional position-switched 
grid-map mode. 

We detected and mapped CO emission in 20 of the galaxies with uniform sensitivity. Here we 
present their global CO properties. The CO emission is generally well centered on the stellar 
disk. However, some Virgo spirals show extended CO distribution to one side or bar-like CO 
distribu tion. The compar ison of our CO data with those of the FCRAO extragalactic CO 
survey (lYoung et al.lll995l ) suggests that the major axis scan and global modeling used by 
Young et al. was largely successful in estimating the total CO luminosity. Young et al. has 
larger CO luminosities than our results for 7 (4) out of 18 galaxies by Icr [2a) uncertainty. 
Further analysis of these data will be presented in our future papers. Extensive multi- 
wav elength data are n ow avail able for joint-analysis with ou r new CO data, including the 
HI fI Chung et allboOTh. radio JYuu. Reddv fc Cond"onll200lh . IR (Spitzer & Akari), optical 
(HST & SDSS), Ha jKoopmann fc Kennevll2004h . and UV (GALEX). We will examine ra- 
dial dependence of the star formation rate and star formation efficiency. Their dependence 
on the morphological type, luminosity, and environments will also be studied using the global 
CO properties. To study environmental effects in different cluster evolutionary stages, we 
will use our Virgo data with the OTF mapping data of Ursa Major cluster spiral galaxies 
(Chung et al. 2010) and single dish spectra of the Pisces filament spiral galaxies (Lee et al. 
2010) obtained using the Kitt Peak 12-m telescope. We will also explore the utility of the 
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CO TuUy-Fisher relation. 

We thank A. Chung for valuable comments and discussions. We would also hke to 
thank T. Jung and H. Kang for excellent support at the data reduction process and the 
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